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Abstract 

In the present work a two-component dark matter model is studied adopting the de- 
generate scenario in the R-parity conserving NMSSM. The gravitino LSP and the neu- 
tralino NLSP are extremely degenerate in mass, avoiding the BBN bounds and obtaining 
a high reheating temperature for thermal leptogenesis. In this model both gravitino (abso- 
lutely stable) and neutralino (quasi-stable) contribute to dark matter, and direct detection 
searches for neutralino are discussed. Points that survive all the constraints correspond to 
a singlino-like neutralino. 



1 Introduction 



There is accumulated evidence both from astrophysics and cosmology that about 1/4 of the 
energy budget of the universe consists of so called dark matter, namely a component which is 
non-relativistic and neither feels the electromagnetic nor the strong interaction. For a review 
on dark matter see e.g. [lj. Although the list of possible dark matter candidates is long (for a 
nice list see [2]), it is fair to say that the most popular dark matter candidate is the lightest 
supersymmetric particle (LSP) in supersymmetric models with R-parity conservation [3]. For 
supersymmetry and supergravity see [I] . The simplest supersymmetric extension of the standard 
model that solves the [i problem [5] is the next-to-minimal supersymmetric standard model 
(NMSSM) [BJ. If we do not consider the axion [7] and the axino [8], the superpartners that have 
the right properties for playing the role of cold dark matter in the universe are the gravitino 
and the lightest neutralino. By far the most discussed case in the literature is the case of the 
neutralino (see the classic review [9]), probably because of the prospects of detection. However, 
in the case in which neutralino is assumed to be the only dark matter component, one has to 
face the fine-tuning problem and the gravitino problem [10] . In most of the parameter space 
the neutralino relic density turns out to be either too small or too large Furthermore, 
unstable gravitinos will undergo late-time cascade decays to a neutralino LSP. These decays will 
destroy the light element abundances built up in BBN, unless Tr < 10 5 GeV [12], which posses 
serious difficulties to the thermal leptogenesis scenario [13]. If, on the other hand, gravitino is 
the LSP and therefore stable, playing the role of cold dark matter in the universe, it is then 
the neutralino that will undergo late time decays into gravitino and hadrons, and the gravitino 
problem is re- introduced [T4"] . 

It has been shown that in the degenerate scenario [15] the BBN and CMB constraints are 
avoided, and high values of the reheating temperature are obtained compatible with thermal 
leptogenesis. Here we focus on the scenario in which the masses of the gravitino LSP and 
neutralino NLSP are extremely degenerate in mass. Under this assumption neutralino becomes 
quasi-stable taking part of cold dark matter of the universe together with gravitino, today is still 
around and can be seen in direct detection searches experiments. 

This article is organized as follows. In the next section we present the theoretical framework. 
In section 3 we discuss all the relevant constraints from colliders and from cosmology, and we 
show our results. Finally, we conclude. 

2 Theoretical framework 

In what folows we review in short the particle physics model, namely the cNMSSM, as well as 
the gravitino production mechanisms. 

2.1 Basics of cNMSSM 

The most straightforward extension of standard model (SM) of particle physics based on SUSY 
is the minimal supersymmetric standard model (MSSM) [IB] . It is a supersymmetric gauge 
theory based on the SM gauge group with the usual representations (singlets, doublets, triplets) 
and on M = 1 SUSY. Excluding gravity, the massless representations of the SUSY algebra are 
a chiral and a vector supermultiplet. The gauge bosons and the gauginos are members of the 
vector supermultiplet, while the matter fields (quarks, leptons, Higgs) and their superpartners 
are members of the chiral supermultiplet. The Higgs sector in the MSSM is enhanced compared 
to the SM case. There are now two Higgs doublets, H u , H^, (or Hi, H 2 ) for anomaly cancelation 
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requirements and for giving masses to both up and down quarks. After electroweak symmetry 
breaking we are left with five physical Higgs bosons, two charged if ± and three neutral A, H, h 
(h being the lightest). Since we have not seen any superpartners yet, SUSY has to be broken. 
In MSSM, SUSY is softly broken by adding to the Lagrangian terms of the form 

• Mass terms for the gauginos g^, Mi, M2, M3 

Mgg (1) 

• Mass terms for sfermions / 

m}Pf (2) 

• Masses and bilinear terms for the Higgs bosons H u , 

m 2 Hu HlH u + m 2 H H\H d + Bfi{H u H d + h.c.) (3) 

• Trilinear couplings between sfermions and Higgs bosons 

AYfiHf 2 (4) 

In the unconstrained MSSM there is a huge number of unknown parameters [17] and thus little 
predictive power. However, motivated by the grand unification idea, the constrained MSSM 
(CMSSM) assumes that gaugino masses, scalar masses and trilinear couplings have (separately) 
a common, universal value at the GUT scale, like the gauge coupling constants do. CMSSM 
is therefore a framework with a small controllable number of parameters, and thus with much 
more predictive power. In the CMSSM there are four parameters, m ,mi/2, A ,tan(3, which are 
explained below, plus the sign of the /i parameter from the Higgs sector. The magnitude of \i, 
as well as the B parameter mentioned above, are determined by the requirement for a proper 
electroweak symmetry breaking. However, the sign of fi remains undetermined. The other four 
parameters of the CMSSM are related by 

• Universal gaugino masses 

Mi(M GUT ) = M 2 (M GUT ) = M 3 (M GUT ) = mifi (5) 

• Universal scalar masses 

mf.(M GUT )=m (6) 

• Universal trilinear couplings 

A% (M GUT ) = A% (M GUT ) = 4 (Mqut) = AoSy (7) 

tan/3 = — (8) 
v-i 

where t>i,t>2 are the vevs of the Higgs doublets and M GUT ~ 10 16 GeV is the Grand 
Unification scale. 
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Unfortunately, the CMSSM suffers from the so-called \x problem [5]. This problem is elegantly 
solved in the framework of the next-to-minimal supersymmetric standard model (NMSSM) [6]. 
In addition to the MSSM Yukawa couplings for quarks and leptons, the NMSSM superpotential 
contains two additional terms involving the Higgs doublet superfields, Hi and H 2 , and the new 
superfield S, a singlet under the SM gauge group SU(3) C x SU{2) L x U(l)y [IE] 



W = tij (Y u H 3 2 Q l u + Y d H{ Q 3 d + Y e Hi V e) - t tj X S H\H{ + \kS z (9) 

where we take Hj = (if", iff ), H 2 = (H^H®), i,j are SU(2) indices, and ei 2 = 1. In this 
model, the usual MSSM bilinear /x term is absent from the superpotential, and only dimensionless 
trilinear couplings are present in W. However, when the scalar component of S acquires a VEV, 
an effective interaction \iHiH 2 is generated, with \i = X(S). 
Finally, the soft SUSY breaking terms are given by [T8] 



— £ so ft = ttlq Q* Q + u* u + d* d + L* L + m^e* e 

2 TT* TT 2 TT* U 1 2 



+ rrf Hl HI Hi + m z H2 H* 2 H 2 + m l s S*S 

+ eij (A u Y u H 3 2 Q'u + Aa Y d HI Q 1 d + A e Y e Hi V e + H.c.) 

+ f -eyA A x SH\Hi + 1 -k A K S 3 + H.c. 

- - (M 3 A 3 A 3 + M 2 A 2 A 2 + Mi Ai Ai + H.c.) (10) 

Clearly, the NMSSM is very similar to the MSSM. Despite the similarities between the 
two particle physics models, the Higgs sector as well as the neutralino mass matrix and mass 
eigenstates in the NMSSM are more complicated and richer compared to the corresponding ones 
in the MSSM. 

In particular, in the Higgs sector we have now two CP-odd neutral, and three CP-even 
neutral Higgses. We make the assumption that there is no CP-violation in the Higgs sector at 
tree level, and neglecting loop level effects, the CP-even and CP-odd states do not mix. We are 
not interested in the CP-odd states, while the CP-even Higgs interaction and physical eigenstates 
are related by the transformation 

K = S ab Hl (11) 

where S is the unitary matrix that diagonalises the CP-even symmetric mass matrix, a, b = 1, 2, 3, 
and the physical eigenstates are ordered as m h o < m h o < m h o. 

In the neutralino sector the situation is again more involved, since the fermionic component of 
S mixes with the neutral Higgsinos, giving rise to a fifth neutralino state. In the weak interaction 

basis defined by ^/° T = (^B° = IU 3 = — i\ 3 , H®, H®, S^j , the neutralino mass terms in the 

Lagrangian are [18] 

4l s = ^(^^° + H.c., (12) 

with M^o a 5 x 5 matrix, 

Mi -M z sin 6 W cos /3 M z sin 6 W sin /3 \ 

Mi Mz cos 9w cos /3 — cos sin /3 

-Mz sin 9 W cos fi M z cos 9w cos /? —As — Xv 2 (13) 

Mz sin 6w sin /J — Mz cos6V sin/3 —As — Xvi 

-\v 2 -Awi 2ks J 



3 



The above matrix can be diagonalised by means of a unitary matrix iV 

^Ai^oN" 1 = diag(m^o, m^o, m^o, m-^o, m^o) (14) 

where m^o is the lightest neutralino mass. Under the above assumptions, the lightest neutralino 
can be expressed as the combination 

Xi = N n B° + N 12 W° + N 13 H° + N lA H° 2 + N 15 S (15) 

In the following, neutralinos with Nf± > 0.9, or Nf 5 > 0.9, will be referred to as bino- or 
singlino-like, respectively. 

Similarly to the CMSSM, in the constrained next-to- minimal supersymmetric standard model 
the universality of mo,Ao,mi/2 at the GUT scale is again assumed, with the only parameters 
now being [T9] 

tan/3, m , A ,m 1/2 , \,A k 
and the sign of the /x parameter can be chosen at will. 

We end the discussion on the particle physics model here, by making a final remark regarding 
the differences between the CMSSM and the cNMSSM. In the CMSSM the lightest neutralino 
is mainly a bino in most of the parameter space, and low values of mo are disfavored because 
they lead to charged sleptons that are lighter than the neutralino x? ; while in the cNMSSM the 
lightest neutralino is mainly a singlino in large regions of the parameter space, thanks to which 
the charged LSP problem can be avoided [19]. Furthermore, in the cNMSSM there are more 
mechanisms that contribute to the neutralino relic density [T9] . 



2.2 Gravitino production 

In the usual case (not in the degenerate scenario) where the neutralino decays into a gravitino and 
standard model particles with a lifetime typically in the range (10 4 — 10 8 ) sec, for the gravitino 
abundance we take the relevant production mechanisms into account and impose the cold dark 
matter constraint [2"U] 

0.1097 < Vt cdm h 2 = n 3/2 h 2 < 0.1165 (16) 

At this point it is convenient to define the gravitino yield, Y3/2 = 713/2/s, where n^/ 2 is the 
gravitino number density, s = /i e //(T)^-T 3 is the entropy density for a relativistic thermal 
bath, and h e ff counts the relativistic degrees of freedom. The gravitino abundance ^3/2 in terms 
of the gravitino yield is given by 

Q 3/2h 2 = G \ 0) 3/2 = 2 ?5 x 1q8 Y- i/2 {T,) (17) 

where we have used the values 

T = 2.73AT = 2.35 x 10~ 13 GeV (18) 

h eff (T ) = 3.91 (19) 

p cr /h 2 = 8.1 x 10~ 47 GeV 4 (20) 

The total gravitino yield has two contributions, namely one from the thermal bath, and one from 
the out-of-equillibrium NLSP decay. 

Y3/2 = Y3/2 + *3/2 P ( 21 ) 
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The contribution from the thermal production has been computed in [211 l22l [23] . In [21] the 
gravitino production was computed in leading order in the gauge coupling g^, in [22J the thermal 
rate was computed in leading order in all Standard Model gauge couplings gy, #2, #3, and in [2"5] 
new effects were taken into account, namely: a) gravitino production via gluon — > gluino + 
gravitino and other decays, apart from the previously considered 2 — > 2 gauge scatterings, b) the 
effect of the top Yukawa coupling, and c) a proper treatment of the reheating process. Here we 
shall use the result of [21] since the corrections of [221 E3] do not alter our conclusions. Therefore 
the thermal gravitino production is given by 

Y*g = 0.29 x 1(T 12 (—^— ) | I ) (22) 



^3/2"— VlO 10 GeV 
or, approximately for a light gravitino, <C m g 




- 10 " 13 lira I - ) < 23 > 



m G 



with itLq the gravitino mass and m g the gluino mass. 

The second contribution to the gravitino abundance comes from the decay of the NLSP 

^/ L 2 SP h 2 = Q NLSP h 2 (24) 

' m NLSP 

with m N LSP the mass of the NLSP, and fl NLS ph 2 the abundance the NLSP would have, had 
it not decayed into the gravitino. In the limit where m NLS p —> tuq and r NLS p 3> 10 17 sec the 
scenario looks as if one would have a two-component dark matter with the NLSP contribution 
^Nisph 2 , and a gravitino contribution from thermal production only, Kwf . Therefore, in the 
degenerate scenario with itinlsp — m c t ne WMAP bound becomes 

0.1097 < Vt cdm h 2 = tt NLSP h 2 + ttlf 2 h 2 < 0.1165 (25) 

where from now on the NLSP is the lightest neutralino, \ = NLSP. 



3 Constraints and results 

- Spectrum and collider constraints: We have used NMSSMTools [21], a computer software that 
computes the masses of the Higgses and the superpartners, the couplings, and the relic density 
of the neutralino, for a given set of the free parameters. We have performed a random scan 
in the whole parameter space (with fixed fi > motivated by the muon anomalous magnetic 
moment), and we have selected only those points that satisfy i) theoretical requirements, such as 
neutralino LSP, correct electroweak symmetry breaking, absence of tachyonic masses etc., and 
ii) LEP bounds on the Higgs mass, collider bounds on SUSY particle masses, and experimental 
data from B-physics [25l [26] . For all these good points the lightest neutralino is either a bino or 
a singlino, and contrary to the case where neutralino is the dark matter particle, here we do not 
require that the neutralino relic density falls within the allowed WMAP range. 

- As we have already mentioned, the total dark matter abundance, and not the neutralino 
one, should satisfy the cold dark matter constraint [20] 

0.1097 < tt cdm h 2 = tt x h 2 + njf 2 h 2 < 0.1165 (26) 
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that relates the reheating temperature after inflation to the gravitino mass as follows 



0.11 = A(m G , rrig)T R + Q x h 2 



(27) 



For a given point in the cNMSSM parameter space, the complete spectrum and couplings have 
been computed, and we are left with two more free parameters, namely the gravitino mass 
and the reheating temperature after inflation. The gravitino mass is equal essentially to the 
neutralino mass, and the precise value can be determined if we specify the neutralino lifetime. 
In the discussion to follow we have used a neutralino lifetime r = 10 26 sec, although the results are 
not sensitive to it, and the figures we have produced for different values of the lifetime cannot 
be distinguished. Finally, the reheating temperature after inflation is obtained from the cold 
dark matter constraint. The thermal production contribution cannot be larger than the total 
dark matter abundance, and for this we can already obtain an upper bound on the reheating 
temperature 



T R < 4.1 x 10 



m 



100 GeV 



TeV 



m-r, 



GeV 



TXIq ~ 



(21 



100 GeV, it 



Assuming a gluino mass m g ~ 1 TeV, we can see that for a heavy gravitino, 
is possible to obtain a reheating temperature large enough for thermal leptogenesis. 

- For neutralino NLSP in the degenerate scenario, the only decay mode is x ~ > for which 
the decay width can be computed once the supergravity Largrangian is known [27], and it is 
given bv pl|l28] 



Nu cos 6 W + N i2 sin 6 



1 w\ 



48ttM2 



m 



m 



1 + 3—^ 

A. 



(29) 



where M* is the Planck mass, m x is the neutralino mass, and 8\y is the weak angle. In the 
limit where the mass difference Am = m x — tjiq is much lower than the masses themselves, 
Am <C m x ,rriQ, the neutralino lifetime becomes 



T = 



1.78 x 10 13 sec (GeV\ 3 
\Nn cos 6w + Nu sin 6w\ 2 \ Am / 

hhflt fnr a mnstlv hinn-nfiiitralinn a l 



(30) 



From this formula one can see that for a mostly bino-neutralino a mass difference of 1 MeV is 
already enough to give a neutralino lifetime larger than the age of the universe. 

- Neutralino-Nucleon spin-independent cross-section: LHC is now running and collecting 
data. Although LHC is a powerful machine to look for physics beyond the standard model, it is 
known that other facilities are also needed to offer complementary information towards the direc- 
tion of searching for supersymmetry and identifying dark matter. The gravitino interactions are 
suppressed by the Planck mass, and therefore direct production of gravitinos at colliders and/or 
direct detection prospects seem to be hopeless. On the other hand, for a weakly interacting 
neutralino there are existing as well as future experiments that put experimental limits on the 
nucleon-neutralino cross- sect ion. The spin-independent cross-section is given by 



Ami 



a X-N 



71 



f. 



where m r is the Nucleon-neutralino reduced mass, m r = m^m x j (thn + m x)> an d 



N 



m N 



E C 

q=u,d,s 



+ 



27 



AN) 
JTG 



E 

q=c,b,t 



"3g 

m a 



(31) 



(32) 
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In the above, f^ G = 1 — ^2 q=u d s f Tq , we have taken the following values for the hadronic matrix 
elements [29] : 

f&l = 0.020 ± 0.004 , f!^ d = 0.026 ± 0.005 , f$ = 0.118 ± 0.062 , 

= 0.014 ± 0.003 , f!$ = 0.036 ± 0.008 , f!$ = 0.118 ± 0.062 . (33) 

and a q is the coupling in the effective Lagrangian 

£ e ff = oii x X Qi Qi (34) 

where i — 1, 2 denotes up- and down-type quarks, and the Lagrangian is summed over the three 
quark generations. The coupling a q can be decomposed into two parts, a q = a q + ot|, where 
the first term is the t-channel exchange of a neutral Higgs (Fig. 1), while the second term is 
the s-channel exchange of a squark (Fig. 2). The expressions for a q is terms of the masses and 
couplings of the model can be found in [T8] . 

Our main results are summarized in the figures below. In Fig. 3 and Fig. 4 we show the 
Nucleon-neutralino spin- independent cross section (in cm 2 ) versus neutralino mass and lightest 
Higgs boson (in GeV) respectively. The blue region corresponds to a bino neutralino, while the 
green region corresponds to a singlino neutralino, and the curves are the current experimental 
limits from CDMS [30] • According to our results the bino scenario is already ruled out, while in 
the singlino case the upper region can be probed by future experiments. In Fig. 5 we show the 
reheating temperature after inflation as a function of the neutralino/gravitino mass. The blue 
region corresponds to a bino, the blue points correspond to singlino, and finally the red points 
correspond to singlino with relatively high values of the cross-section, namely <j x -n > 10~ 47 cm 2 . 
The largest values of Tr correspond to a bino, which is ruled out, and for the singlino with 
relatively high values of cross-section we obtain a reheating temperature Tr ~ 5 x 10 9 GeV for 
a neutralino/gravitino mass m x ~ tuq ~ 200 GeV . In the last figure we show the (mo-mi/2) 
plane (mo and m.1/2 in GeV) for singlino points with a cross-section larger than 10 -47 cm 2 , or 
lower than 10~ 47 cm 2 . We see that m is not larger than 600 GeV, and therefore future direct 
detection experiments cannot probe a region of the parameter space which can neither be probed 
by LHC. 

Before ending the discussion, let us briefly comment on a possible collider signature of our 
model. In the singlino-like neutralino case, where the coupling A is small, A < 0.01, the lightest 
Higgs can be very light, m^ < 60 GeV, which has a significant singlet composition, thus escaping 
detection and being in agreement with accelerator data. In this case the next-to-lightest Higgs 
is the SM-like Higgs, with a mass m# ~ (116 — 118) GeV, and the decay channel H — > hh is 
kinematically allowed. Since the lightest Higgs is expected to exit the detector without been 
seen, the decay channel H — > hh is an invisible one. This is to be contrasted with the cases of 
SM and of the MSSM, where the Higgs (in the SM) and the SM-like Higgs (in the MSSM) with 
a mass in the above range decays (almost entirely) into bb and t + t~, with the sum of the two 
branching ratios being practically one jSTJ. In Table 1 we show the range of the parameters of 
the model where we obtain a very light Higgs and the decay channel H — > hh is kinematically 
allowed. For most of the points the branching ratio is negligible, even as low as ~ 10 ~ 9 , but 
points exist for which the branching ratio becomes sizable, BR ~ 0.1, with the maximum value 
obtained being BR(H — > hh) = 0.13. In Table 2 we show the Higgs boson masses Mfj, and 
the branching ratio for four points. 
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A 


tanfi 


A (GeV) 


A k (GeV) 


m (GeV) 


M 1/2 (GeV) 


0.0102 


34.89 


-127.24 


-106.41 


124.28 


606.5 


0.000127 


25.49 


-225.09 


-170.8 


10.6 


401.43 



Table 1: Range of the parameters of the model where the decay channel H — y hh is kinematically 
allowed. The last row shows the minimum value, while the row in the middle shows the maximum 
value of the parameters. 



A 


tan/3 


A (GeV) 


A k (GeV) 


m (GeV) 


M 1/2 (GeV) 


M H (GeV) 


M h (GeV) 


BR(H hh) 


0.000107 


34.69 


-202.66 


-115.23 


22.65 


581.2 


116.97 


54.81 


3.4 x 10~ 9 


0.00216 


30.25 


-158.28 


-135.98 


49.33 


467.87 


115.8 


46.8 


2.82 x lO" 4 


0.009154 


32.3 


-181.5 


-131.01 


16.2 


540.57 


118.33 


49.44 


0.13 


0.005064 


30.94 


-185.46 


-169.65 


35.86 


552.87 


117.29 


57.17 


0.01 



Table 2: Higgs boson masses, the branching ratio and the values of the parameters for four 
points. 



4 Conclusions 

In the framework of NMSSM, which solves the /z problem, we have assumed that the gravitino 
LSP and the lightest neutralino NLSP are degenerate in mass. Under this assumption the 
neutralino becomes extremely long-lived avoiding the BBN bounds. In this scenario we have 
a two component dark matter made out of the absolutely stable gravitino and the quasi-stable 
neutralino. We have performed a random scan over the whole parameter space keeping the points 
that satisfy the available collider constraints plus the WMAP bound for dark matter. These 
points correspond to either a bino or a singlino neutralino. We have computed the neutralino- 
nucleon spin-independent cross section as a function of the neutralino mass and the lightest 
Higgs mass, and we find that the bino case is ruled out (see Fig. 3 and Fig. 4). Then we 
explored the (m — mi/ 2 ) parameter space, and the reheating temperature dependence of the 
neutralino/gravitino mass for the singlino points that correspond to cross section values to be 
probed by future experiments. Finally, we have briefly discussed an interesting possibility for 
collider signatures, namely the possibility of having an invisible decay channel H — y hh, where H 
is the SM-like Higgs and h is the lightest Higgs that escapes detection, with a sizable branching 
ratio and maximum allowed value BR(H — y hh) ~ 0.13. 
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Figure 3: Spin-independent neutralino-nucleon (proton) cross-section versus neutralino mass. 
Shown are the available experimental bounds from CDMS, and the predictions of the theoretical 
model. The blue region corresponds to a bino-like neutralino, while the green points correspond 
to a singlino-like neutralino. 
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Figure 4: Spin-independent neutralino-nucleon (proton) cross-section versus the lightest Higgs 
mass. Shown are the available experimental bounds CDMS, and the predictions of the theoretical 
model. The blue region corresponds to a bino-like neutralino, while the green points correspond 
to a singlino-like neutralino. 
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Figure 5: Reheating temperature versus neutralino (or gravitino mass). Blue points correspond 
to bino, green points correspond to singlino, and red points correspond to singlino with a cross- 
section larger than 10~ 47 cm 2 . 
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$ure 6: The (7710-7711/2) plane for the singlino points. One color corresponds to a cross-section 
ger than 10~ 47 cm 2 , and the other color corresponds to a cross-section lower than 10 



cm . 
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